Exact low-temperature properties of a class of highly frustrated Hubbard models 
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We study the repulsive Hubbard model both analytically and numerically on a family of highly 
frustrated lattices which have one-electron states localized on isolated trapping cells. We construct 
and count exact many-electron ground states for a wide range of electron densities and obtain 
closed-form expressions for the low-temperature thermodynamic quantities which are universal for 
all lattices of the family. Furthermore, we find that saturated ferromagnetism is obtained only for 
sufficiently high electron densities and large Hubbard repulsion U while there is no finite average 
moment in the ground states at lower densities. 
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I. INTRODUCTION 

Exactly solvable interacting quantum lattice sys- 
tems are of great importance in condensed matter 
physics. Although most of the examples are known 
for one-dimensional systems, two-dimensional and three- 
dimensional models are also coming into sight. A cru- 
cial point in searching for exactly solvable lattice models 
concerns the lattice geometry. Furthermore the inter- 
play between lattice geometry, interactions and quantum 
fluctuations often gives rise to exotic quantum phases. 
Famous examples where a special arrangement of in- 
teraction bonds allows to find exact quantum many- 
body ground states are the Majumdar- Ghosh model, ^ the 
Shastry-Sutherland model^, and the Kitaev modeli^ 

Besides revealing new magnetic properties of previ- 
ously known materials (e.g. the diamond-chain com- 
pound azurite^) there are presently various possibilities 
to design interacting lattice systems with controlled ge- 
ometry. Modern strategies in chemistry open a route to 
synthesize new materials with a desired lattice structure 
and intersite interactions J^i^ Moreover, recent progress 
in nanotechnology allows the fabrication of quantum dot 
superlattices and quantum wire systems with any type of 
latticcji Another rapidly developing field is the controlled 
setup of optical lattices for cold atomsi^ii^ 

Motivated by these achievements we propose here a 
class of lattices (including the well-known diamond chain, 
frustrated ladder, square-kagome and checkerboard lat- 
tices) for which various properties of the Hubbard model 
can be examined rigorously. In particular, we character- 
ize the complete manifold of highly degenerate ground 
states for electron numbers n = 1, . . . ,J\f (x N and calcu- 
late low-temperature thermodynamic quantities around 
a particular value of the chemical potential fiQ. The gen- 
eral lattice construction rules illustrated below are based 
on a local point of view similar (but not identical) to ear- 
lier considerations for elcctronic^*''^^ and spin systems^— 
Below we also discuss some properties of the Hubbard 
model for a few one-dimensional and two-dimensional 
representatives. We illustrate the dominating role of the 



exact ground states for the low-temperature physics of 
the Hubbard model at certain electron densities. Fur- 
thermore we analyze the ground state with respect to 
magnetic properties. We confirm our analytical findings 
by numerical data for finite systems. 



II. LATTICES WITH TRAPPING CELLS 



We consider the iV-site Hubbard Hamiltonian 



H= Y, Hoa+Hu, Hu^uY,' 



'■iA^idy 



^=^1 



Hoa = Y ^^'J (4o-Cj\<T + cl^C^^a) + M XI '*'^'^' (^) 



(«J> 



where i denotes the lattice sites, {i,j) denote the bonds 
connecting neighboring sites, the cj^ {ci,cr) are the usual 



fermion operators. 



J 



ti,j > are the hop- 



ping integrals, [/ > is the on-site Coulomb repulsion, 
and /i is the chemical potential. 

We consider the Hubbard model ([T|) on a family of 
lattices defined by the following construction rules: (i) 
Take a 'trapping cell', i.e., a finite region where an elec- 
tron in the infinite lattice will be localized. For sim- 
plicity, we consider here bipartite cells with equivalent 
sites and bonds, such as a single bond between two sites, 
equilateral even polygons, or a cubet^ (ii) Solve the 
one-electron problem for the trapping cell, finding the 
lowest-energy eigenfunction oc J^i'^i'^l a\^) with |0) de- 
noting the vacuum state. For traps consisting of a single 
bond or a square the ground state is nondegenerate and 
ai = —02 = 1 or fli = —02 = a3 = -04 = 1, respectively. 
(iii) Arrange the trapping cells into a regular pattern so 
that trapping cells do not touch each other (do not have 
common sites) and complete the lattice by connecting the 
cells via surrounding (connecting) bonds. Most impor- 
tantly, the connecting bond scheme should prevent the 
escape of the localized electron from the trap, i.e., the 
constructed one-electron (localized) state should remain 




FIG. 1; (Color online) Diamond chain (a), frustrated ladder 
(b), square- kagome (c), and checkerboard (d) lattices. The 
circles indicate the lattice sites, the thick (thin) lines indicate 
the hopping paths along the trap (connecting) bonds. For 
the ideal diamond chain with fi — ta we will use the notation 
t — ti = t^, t — t2. 



an eigenstate of the Haniiltonian on the infinite lattice. 
It is easy to show that a sufficient condition for this is 
X]j tr,iCii = 0, where the sum runs over all sites z of a 
trapping cell and r is an arbitrary site which does not 
belong to the trap, see also Ref. [l^. For such traps as 
a single bond, square or cube this condition is fulfilled if 
an arbitrary bond of the trapping cell and the surround- 
ing bonds attached to the two sites of this bond form an 
isosceles triangle, (iv) Choose the hopping integrals t' of 
the trapping cells sufficiently large, t' > t'^, so that the 
localized states become the lowest-energy ones in com- 
parison with other (extended) one-electron states. 

Following the rules formulated above we are able 
to construct many different one-dimensional, two- 
dimensional, and even three-dimensional lattices. Note 
that a lattice constructed in this manner has a flat one- 
electron band which for t' > t'^ becomes the lowest-energy 
one. However, we emphasize that such lattices do not be- 
long to Tasaki's models for flat-band ferromagnetism.^" 

In what follows we focus for concreteness on some typ- 
ical representatives which were investigated in the con- 
text of strongly correlated systems earlier, namely, the 
diamond chain^i^ the frustrated ladder ;^^'^^ the square- 
kagome, and checkerboard lattice o^^'^^ (see Fig. [1]). It 
is also convenient to denote sites by to,p, where m — 
1, . . . , AA denotes unit cells {N = A^/3, N/2, N/&, N/A 
for the diamond chain, frustrated ladder, square- kagome, 
and checkerboard lattice, respectively) and p denotes 
sites inside the unit cell (see Fig.[T]). Note that for the lat- 
tices considered here the number of unit cells is identical 



to the number of trapping cells. However, the trapping 
cell may or may not be identical to the unit cell. 



III. LOCALIZED-ELECTRON STATES FOR 
[7 = AND U >Q 

Next we turn to the characterization of the complete 
manifold of highly degenerate ground states for electron 
numbers n = 1,...,7V oc TV. We start with the case 
[/ = considering the diamond chain with t ^ ti = t^, 
t' = t2 for concreteness. Three one-electron bands. 



El- fl 



-t', 



£2.3(K)-/i = |ty(|) +4t2(l + cosK), (2) 

are arranged as follows for t' > t'^ = 2t: £i < £2('«) < 
e3{K). Note that the lowest one-electron band ei is 
completely flat, i.e., K-independent. The correspond- 
ing states can be localized in real space. We introduce 



the operators ij^ 



^m,2,a 






[indices 2 and 3 



denote the bottom and top sites on the vertical bond 
(see Fig. [Ij, m enumerates the unit cells] which satisfy 
[TJocTj^m.CT]- — £i ^m cr- Then all 2A/' one-electron states 
belonging to the flat band can be written as lj^^\Q). 
Application of n distinct operators l}^ ^ to |0) yields n- 
electron states with energy En = nsi, n = l,...,Af. 
Note that all trapping cells are disconnected and the de- 
generacy of these n-electron states is 



,(0) 



9):}'in) 



f2N 



(3) 



These arguments can be applied to other models: for the 
energy of the flat band ei — /i we find — t', — 2t', — 2t' for 
the frustrated ladder, square-kagome, and checkerboard 
lattice, respectively. We assume t' /t > t'^/t = 2, 1, 1 for 
the frustrated ladder, square-kagome, and checkerboard 
lattice, respectively. For the square-kagome and checker- 
board lattice an electron may be localized on smallest- 
area squares and Z„ ,j = c^^^^ — c„^2.a + ^m,?,,^ ~ ^m,4.a^ 
where the indices 1, . . . , 4 denote the vertices of the 
square. 

We now address the case U > 0. Since Hu is a positive 
semidefinite operator for L'^ > 0, it can only increase ener- 
gies. The states for which each trapping cell contains up 
to one electron are exact eigenstates of the full Hamil- 
tonian ([T]) with the [/-independent energy En = nei 
and thus they remain the ground states in the subspaces 
n = 2, ... ,M in the presence of a Hubbard repulsion 
U > Q. It is straightforward to count these localized 
n-electron ground states and we find 



gMn)^r'(^)<9]^\n). 



(4) 



The localized n-electron states are linearly indepen- 
dent which can be proven using the arguments of Ref. Il9| 



(orthogonal class in the nomenclature of that reference) . 
These states are the only ground states in each subspace 
n — 1 , . . . , A/". This can be seen by recalling from spin 
systems that a finite separation of the flat one-particle 
band from the next (dispersive) band ensures complete- 
ness of the localized states.— For the present models, we 
can even control the energy gap by varying i'/t, in con- 
trast to the sawtooth Hubbard chain.- - The state with TV 
electrons possesses perfect charge ordering (each trapping 
cell is occupied by precisely one electron) and therefore it 
can be understood as a particular realization of a Wigner 
crystal.^ However, with respect to the spin orientations 
it has a huge degeneracy 2^ (compare the degeneracy of 
the ground state, J\f + 1, at quarter filling when n = J\f 
for the sawtooth chai n^^i^'^ ). 



C{T, /i, N). Eqs. (5) - (8) given above represent the con- 
tribution of the localized-electron states to the respective 
thermodynamic quantities. In section V we will demon- 
strate that these analytic expressions are an excellent de- 
scription of the thermodynamics of the full model at low 
temperatures. 

The average ground-state electron number n(0, fj-o, N) 
exhibits a jump from TV to zero as fi exceeds fiQ. More- 
over, at /i — /io we have x = resulting in a finite resid- 
ual entropy s[o, fj.o,N)/J\f = ln3 « 1.0986. We note 
that thermodynamic quantities for all considered models 
are identical up to a factor M/N and the concrete value 
of Ho, i.e., the low-temperature behavior is universal for 
the whole family of lattices constructed by the rules for- 
mulated in section im 



IV. THERMODYNAMICS 



V. NUMERICAL RESULTS 



The localized-electron ground states in the sectors 
with n < J\f have important implications for the low- 
temperature properties of the Hubbard model ([T]) around 
a chemical potential /io = ^^ — £i'- Due to their huge 
degeneracy they dominate the grand-canonical partition 
function at low temperatures. Knowing the degeneracy 
dH) of the ground states guin) and their energy En = nei 
we can write the grand-canonical partition function as 



A/" 

5(r, M, TV) = ^ g^in) e^^^/^ = (l + 2 e"^^/^ 



n=0 



H 



(5) 

Note that only the combination x = —ei/T = {fio — ^j.)/T 
enters the localized-electron contribution to any ther- 
modynamic quantity. The thermodynamic potential be- 
comes n{T,^i,N)/Af ^ -Tln(l + 2expa;) (the rhs. is 
valid for both finite TV and TV -^ oo) leading to simple 
expressions for thermodynamic quantities such as 






2e^ 



1-H2e^ 



(6) 



for the average electron density n{T,fj.,N) 

dn{T,ti,N)/dfi, 

^(^''^'^)^ln(l + 2e^) '^^^ 



M 



l-|-2e^ 



(7) 



for the entropy S{T, fj,, TV) = -(9fi(r, fj,, N)/dT, or the 
specific heat in the grand-canonical ensemble 



C(T,/i,TV) _ T^S{T,^i,N) 



N 



MdT 



2x^e^ 
l + 2e^ 



2.Te^ 
l + 2e^ 



(8) 



The specific heat at constant n vanishes identically, 
C{T,n,N) = 0, as can be verified explicitly using the 
results ©, dZD, and ^ for n{T,^i,N), S{T,fj,,N), and 



Let us now present numerical results obtained by ex- 
act diagonalization for finite lattices. We set t = 1 for 
convenience. These calculations have been performed, 
on the one hand, to estimate the range of validity of the 
derived expressions for the low-temperature thermody- 
namics and, on the other hand, to study the states with 
n > TV in more detail. The numerical effort to diagonalize 
the Hubbard model H]) grows rapidly with TV. It is there- 
fore convenient to consider also the limit [T — > cxd where 
doubly occupied sites can be eliminated from the Hilbert 
space. Computations were performed using different pro- 
grams including J. Schulenburg's spinpack^*' and a cus- 
tom implementation of the Householder algorithm. ^-'^ 

Note first that the degeneracies of the localized ground 
states calculated for various lattices and parameter sets 
perfectly fit to the prediction gj^f{n) of section Hill The 
average electron density n/TV versus /i//io shown in Fig. [2] 
exhibits a jump between and 1 at /i//io = 1 and a 
plateau at n/Af = 1 with the width A/z//i. The charge 
gap A^ = E{M+ 1) - 2E{M)+E{M-l) determines the 
region in which the localized states exclusively control 
the ground-state behavior of the model ([T]). The inset 
in Fig. [2] shows that A/i is almost size- independent. The 
charge gap increases almost linearly with U for small [/, 
showing that this is a correlation effect. For the diamond 
chain with t' — 3t, the charge gap saturates at A/i k, 
0.43 1 for [T > 3 1. 



A. Thermodynamics for ideal geometry 

Next we compare the specific heat C{T,ii,N) in the 
grand-canonical ensemble for the diamond chain, the 
frustrated ladder, the square-kagome lattice, and the 
checkerboard lattice at /i = 0.95 /io, Mo, and 1.05 /io 
in the limit U ^ oo (Fig. [3]). C{T,n,N)/M versus 
r//io exhibits a universal additional low-temperature 
maximum for the chemical potential around /^o- This 
low-temperature maximum emerges due to the manifold 
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FIG. 2: (Color online) Electron density njM versus chemical 
potential /i//io for the diamond chain (A'' = 18, i' = 3), the 
frustrated ladder (A'^ = 16, t' = 3), the square-kagome lattice 
(iV = 24, f' = 2), and the checkerboard lattice (iV = 16, 
t' = 2) for [/ — > cx) and T = 0. The universal dependence (|6]) 
at r = is given by 9{1 — fi/fio)- Inset: Charge gap Afi/fio at 
n/A/" = 1 versus U for the diamond chain [A^ = 12 (triangles), 
A- = 18 (line), i' = 3]. 
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of localized-electron states, whose energies are slightly 
split for /i = 0.95 /io and 1.05 /iq but still well sepa- 
rated from energies of other higher energy states. In- 
deed, the low-temperature maximum is excellently de- 
scribed by the localized-electron formula ([5]) (lines in 
Fig. [S]). Recall that the result dH) is valid for any N. 
Accordingly, there are no finite-size effects for the low- 
temperature maximum of the specific heat. In Fig. [31 
one observes another maximum of C{T, fi, N) at higher 
temperatures. This maximum collects all states which 
are not localized-electron states and thus depends on de- 
tails of the model. If /.t — > /ig the low-temperature maxi- 
mum shifts to lower T and is even much better separated 
from the high-temperature maximum. At fi ~ fj,o, we 
have X — 0. Consequently, the contribution ([5]) of the 
localized-electron states to the specific heat C vanishes 
identically. Therefore, the temperature where the numer- 
ical data for C{T, fio,N) begin to deviate from zero (see 
the middle panel of Fig. [3]) represents a characteristic 
temperature below which the localized states exclusively 
control the thermodynamic behavior of the model H]). 

Fig. m illustrates the influence of different values of 
U, using the example of the frustrated ladder with 
A'^ = 12 and t' — 3. The low-temperature maximum 
in C{T, /io, N) is independent of U, as expected. By con- 
trast, the high-temperature maximum does depend on 
U and it extends to lower temperatures for smaller val- 
ues of U. Nevertheless, the localized-electron states are 
still well separated from the other states for the values of 
[/ > 4 and fi = 0.95 /xq, 1-05 /iq shown in Fig. H) 

Fig. [5] shows the entropy S for the same lattices and 
parameters as in Fig. [31 Since C{T, n, N) is a temper- 
ature derivative of S{T, ii,N), the features in the latter 
quantity correspond to those in the former. In partic- 



FIG. 3: (Color online) C{T, iJ.,N)/J\f versus temperature 
T/fio in the limit [/ — > oo for the diamond chain (triangles), 
the frustrated ladder (squares) , the square-kagome lattice (di- 
amonds), and the checkerboard lattice (circles). We also show 
the universal dependence (jS)! (lines). 
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FIG. 4: (Color online) C{T, fj., N)/M for the frustrated ladder 
with A' = 12 and t' = 3 for different values of U . We also 
show the universal dependence ([H (lines). 
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FIG. 5: (Color online) Entropy S{T,^,N)/M versus temper- 
ature r/^o in the limit [/ — > cx3 for the diamond chain (tri- 
angles), the frustrated ladder (squares), the square-kagome 
lattice (diamonds), and the checkerboard lattice (circles). We 
also show the universal dependence (O (hues). 



iilar a steep slope in S{T, fi, N) corresponds to a max- 
imum of C{T, fi,N). However, there is one additional 
piece of information in S(T^ fi, N), namely its value in 
the low-temperature limit: for /i = /iq, S{T, fi, N)/J\f 
tends to In 3 « 1.0986 for T ^ (see section |IV1) . If the 
chemical potential /i is in the charge-gap region, i.e., for 
Ho — A/i < /J, < /io, the entropy per unit cell stays also 
finite for T ^ (compare the upper panel of Fig. [5]) . For 
these values of /x the ground state is the charge-ordered 
state (Wigner crystal) with n = M electrons and accord- 
ingly exhibits a ground-state degeneracy oi gj\f{M) = 2^ 
due to the independence of the spin orientations of the 
J\f electrons, see Eq. (U)). Indeed, Eq. ([7]) shows that the 
entropy per cell S/N — > ln2 w 0.6931 when x — > oo, i.e., 
T ^ and /i < /io- 



B. Deviations from ideal geometry 

In real compounds the conditions under which the lo- 
calized states exist may be not strictly fulfilled. It is 
therefore important to study deviations from 'ideal ge- 
ometry'. Such deviations or distortions in general lift the 
ground-state degeneracy. Nevertheless, one may expect 
that localized-state effects survive in a certain tempera- 
ture range if the distortion is sufficiently small such that 
the originally degenerate energy levels remain close to 
each other. We investigate this issue in more detail using 



FIG. 6; (Color online) C{T,fi,N)/J\f for the ideal diamond 
chain {t' = 4, f = 1) and the distorted diamond chain (ii = 
1.1, f2 = 4, fa = 0.9). The system size is A'^ = 12 in all cases. 



the example of the distorted diamond chain {N — 12, 
t' — 4) where we distinguish the hopping terms running 
from north-west to south-east ti > t (ti = 1.1) and those 
running from south-west to north-east t^ < t {t^ — 0.9), 
see Fig. [1] This generalization is inspired by the set of ex- 
change interactions proposed for azurite in Ref. 0- Fig. [6] 
shows numerical results for N = 12, U = oo and 10 in 
comparison to the data for the undistorted case which 
are identical to those shown in Fig. [31 The effect of the 
distortion is evidently very small for T > 0.04. In par- 
ticular, the additional low-temperature maximum in the 
grand-canonical specific heat C{T,fi,N)/M for ^ = 3.8 
and 4.2 (corresponding to ^ = 0.95 /iq and 1.05 /iq in the 
undistorted case) is essentially unaffected. This implies 
that an important fingerprint of the highly degenerate 
localized states survives a small distortion. Additional 
features emerge in Fig. [S] for the distorted diamond chain 
in the region T < 0.01 at /i = 3.8 and for T < 0.1 
at /x = 4. It is evident from Fig. [S] that these features 
depend on the value of \J . Comparison with smaller sys- 
tem sizes for [/ = 10 (not shown) exhibits also finite- 
size effects for T < 0.001 at ^ = 3.8 and for T < 0.01 
at ^ = 4. One may therefore speculate that the addi- 
tional low-temperature maximum observed in C{T, fj,, N) 
for // = 4 around T « 0.006 in Fig. [6] survives the ther- 
modynamic limit, but it is difficult to infer the behavior 
in the thermodynamic limit at even lower temperatures. 
The specific heat C{T, n, N) in the canonical ensem- 
ble, i.e., for a fixed number of electrons n, may be par- 
ticularly relevant from the experimental point of view. 
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FIG. 7: (Color online) Specific heat C{T,n,N)/Af for the 
ideal diamond chain {t' — 4, t — 1) and the distorted diamond 
chain (ii = 1.1, ^2 = 4, ta = 0.9) with n = N/3 electrons for 
iV = 9, 12, 15, 18. The Hubbard repulsion is [/ = 10. 



In section IIVI we have already used the results for the 
grand-canonical ensemble to point out that the contribu- 
tion of the localized-electron states to C{T, n, N) always 
vanishes identically in the ideal situation. This is also ev- 
ident if one considers a fixed number of electrons: in this 
case all localized n-electron states have the same energy 
E and using the representation C = {{E"^) — {E)"^) /T^ 
in terms of the fiuctuations of E one finds C{T, n, N) = 0. 
We use again the example of the (distorted) diamond 
chain to take a closer look at the behavior of C{T, n, N). 
The case of the Wigner crystal, i.e., n = M may be 
particularly interesting. Fig. [7] shows numerical results 
for the specific heat C(r, n, N)/Af of both the ideal and 
the distorted diamond chain with n ~ Af — N/3 elec- 
trons. First, we observe that C{T,n,N) of the ideal 
diamond chain (solid lines in Fig. [7]) is indeed indistin- 
guishable from zero for T < 0.1 and accordingly there 
is no low-temperature maximum. By contrast, for the 
distorted diamond chain the exact degeneracy of the 
ground-state manifold is removed, leading to the recov- 
ery of a low-temperature maximum. For the parame- 
ters used in Fig. [7l this low-temperature maximum ap- 
pears at T w 10"'' (see dashed lines in Fig. [7|), reflecting 
the fact that the detuning ^3 = 0.9 7^ ti = 1.1 leads 
only to a very small splitting of the originally 2"-fold de- 
generate ground-state manifold. There is another high- 
temperature maximum at temperatures of order one in 
Fig. [7] where the results for the distorted diamond chain 
are indistinguishable from those for the ideal situation. 
The double-peak structure found in the canonical specific 
heat C{T, n, N) of the distorted situation is qualitatively 
remarkably similar to the double-peak structure in the 
grand-canonical specific heat C{T, ^, N) observed e.g. for 



^ = 3.8 in Fig. H 

While finite-size effects are mostly unimportant in the 
grand-canonical ensemble, they are clearly relevant in 
the low-temperature region of Fig. [7] and visible even 
in the region of the high-temperature maximum. Such 
enhanced finite-size effects are typical for computations 
performed in the canonical ensemble where one considers 
only a fixed number of electrons n. 



VI. MAGNETIC PROPERTIES 

We return now to the ideal geometry and discuss 
ground-state magnetism. Since the trapping cells are in- 
dependent of each other (i.e., they do not have common 
sites), the cells can be occupied independently by elec- 
trons with spin up or spin down for n = 1, ... ,7V. As a 
result, for n < Af the majority of the degenerate ground 
states of the Hamiltonian ^ on the considered lattices 
is nonmagnetic. Indeed, using the localized-electron pic- 
ture of section lllli it is straightforward to show for the 
total spin operator S = (S"^, S"^, S^) that 



iV2 iV2 4iV2 ' 







(9) 



for n < Af , where the limit A^ -^- 00 is taken for a fixed 
electron density n/N . Here (. . .)„ denotes the average 
over all ground states in the sector with n electrons. 

The situation changes fundamentally aX n = Af +1. If 
U is small (with respect to the energy gap between the 
flat band and the next band; this gap can be controlled 
by the ratio t' /t) an extra electron may lead to formation 
of a complicated many-body state with an energy that in- 
creases with growing U . However, if U is large enough, 
U > Uc[Af + 1), it might be energetically favorable to 
avoid double occupancy by putting one extra electron 
in the lowest-energy state of the next (dispersive) one- 
electron band in addition to the Af localized electrons. 
The magnitude of Uc depends on the system under con- 
sideration. It is possible to find an explicit expression for 
such an eigenstate: 



?t ,t 



Wn+i) «/3lo,T^AA,T---h,T|0> 



(10) 



Here (3\^ „ creates an electron in the lowest dispersive 
band with momentum kq and spin a. The state \(pj\f+i) 
is fully polarized, i.e.. 



{Lpj^+l\S'^\ipj^+l) 



N+l fAf + l 



(11) 



and it has a t/-independent energy Afei + £2(^0). Other 
states belonging to a spin- ^""^ SU(2)-multiplet can 

be obtained by applying S~ — X^iS |Ci,T to this state 
(Kramers degeneracy). This multiplet can also be ob- 
tained by applying 5+ = ^^ c| ^Ci^i to the spin down 
counterpart of Eq. (jTU]) . 



TABLE I; Number of electrons n > 1 for which there is 
a unique saturated ferromagnetic ground-state multiplet in 
the limit [/ — > cxa on the specified lattice subject to peri- 
odic boundary conditions. Note that saturated ground-state 
ferromagnetism for the considered lattices of N sites occurs 
only for a sufficiently large number of electrons J\f < n < N 
(apart from the trivial case n = 1). Note that for the 
square-kagome lattice with N = 24 sites the ground states 
for n = 16, 17, 18, 19 have not been calculated, because of the 
large size of the Hamiltonian matrix. 



lattice 


number of electrons 


diamond, iV = 12 (A^ = 4) 


n = 5,7,8,9,11 


diamond, N = 18 {Af = 6) 


n = 7,9,11,12,13,15,17 


ladder, iV = 12 (A^ = 6) 


n = 7,9, 11 


ladder, iV = 16 (A^ = 8) 


n = 9,11,13,15 


square-kagome. A*' = 18 (TV = 3) 


n = 4,5,7,8, 15, 17 


square-kagome, A = 24 {Af = 4) 


n = 7, 8, 9,..., 20, 21, 23 


checkerboard, A = 16 (A^ = 4) 


n = 7,8,9, 15 



Numerically, we find such ferromagnetic ground states 
for large U also for bigger M < n < N. For finite sys- 
tems the numbers of electrons n for which ferromagnetic 
ground states appear depend on the lattice size and the 
boundary conditions i^^'^'^ For all lattices considered in 
this paper and the imposed periodic boundary conditions 
we find fully polarized ground states for particular val- 
ues of the electron number n > J\f for sufficiently large 
U > Uc- In Table [l] we list some combinations of finite 
lattices and numbers of electrons n for which the ground 
state is fully polarized for U ^ oo. Note that the state 
pop is incompatible with the boundary conditions of the 
A^ = 24 square-kagome and the iV = 16 checkerboard 
lattice such that there is no saturated ground-state fer- 
romagnetism for n = 5 on these two lattices. In addition, 
we give as an example the strength of correlation U > Uc 
that is needed to realize ground-state ferromagnetism for 
n = M + 1 for the frustrated ladder with iV = 16 sites: 
we found Uc « 0.42 and 12.78 for t' = 2.1 and 4, respec- 
tively. 

This kind of ferromagnetism that occurs only for suf- 
ficiently strong on-site repulsion (note that it is not 
fiat-band ferromagnetism that occurs for any U > 
0) was discussed also earlier for some one-dimensional 
systems jiii^Sii^^ Moreover we mention that the satu- 
rated ferromagnetism found for the largest electron num- 
ber n = N — 1 corresponds to Nagaoka's well-known 
result F^^i^i However, we emphasize that the occurrence of 
ferromagnetism for other electron numbers in the range 
JV < n < N is a generic feature of the considered lat- 
tices and is therefore not restricted to one-dimensional 
systems. 



VII. CONCLUSIONS 

In summary, we have given an exact solution for the 
ground-state properties of a correlated many-electron 
system on a class of lattices in a certain range of 
the chemical potential /i. We have studied the low- 
temperature thermodynamics which for fi around /iq is 
controlled just by the manifold of localized ground states. 
In particular, we have presented explicit expressions for 
the low-temperature behavior of the grand-canonical par- 
tition function and related quantities such as the grand- 
canonical specific heat C{T, fi, N). The localized-electron 
features have no finite-size effects and are universal, i.e., 
they are the same for the whole class of lattices. 

Apart from the high-temperature maximum that is 
typical for systems with a bounded energy spectrum 
we find an additional low-temperature maximum in 
C{T,fi,N) if the chemical potential /i deviates slightly 
from /iQ. This extra maximum survives under small de- 
viations from ideal lattice geometry. Moreover, for the 
canonical specific heat C{T, n, N) the lifting of degener- 
acy caused by deviations from ideal lattice geometry ac- 
tually gives rise to an additional low-temperature max- 
imum in the specific heat. Thus, any splitting of the 
highly degenerate localized electron states, be it by a 
deviation of the chemical potential from hq or by a devi- 
ation from ideal lattice geometry, leads to an extra low- 
temperature maximum in the specific heat as a charac- 
teristic fingerprint of the localized-electron states. 

Furthermore, we argued that in contrast to the flat- 
band ferromagnetaiS there is no ground-state ferromag- 
netism in the present class of models if the electron den- 
sity satisfies n/JV < 1. For electron numbers M < n < 
N we observe saturated ground-state ferromagnetism 
for all considered lattices, including square-kagome and 
checkerboard lattices. An explicit expression for the sat- 
urated ferromagnetic ground state in the sector n = A/'+l 
is given in Eq. (fTO|) . 

Finally we note that the t— J model on the considered 
lattices may exhibit similar localized ground states in the 
subspaces with n = 1, . . . ,M electrons for values of the 
spin exchange interaction J up to about the hopping in- 
tegral t. Another interesting variation on the considered 
models refers to multiorbital systems.^ For instance, the 
one-orbital Hubbard model on a frustrated ladder can be 
related to a two-orbital Hubbard model on a simple chain 
with a hybridization term corresponding to the hopping 
on a rung. This relation between one-orbital models and 
multiorbital models extends the region of applicability of 
the localized-states scenario. 
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